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In momentum space the Feynman propagator Dp{k) at non-zero temperature is de- 
fined by a simple dispersion relation with the familiar property of being an even function 
of and analytic for Re(/c^)^ > 0. The coordinate space form of the propagator Dp{x) 
is expressed directly in terms of matrix elements of the field operator and requires a new 
type of operator ordering. 
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I. INTRODUCTION 

In zero-temperature field theory the Feynman propagator is the vacuum expectation 
value of the time-ordered product of two field operators. Therefore it is rather curious 
that at non-zero temperature the Feynman propagator and the time-ordered propagator 
propagator are different. At T the time-ordered propagator for a real scalar field is 

„— /3-Ejv 

iD,,ix) = Y,{N\T{Hx)m)\N) (1) 

AT 

where |A^) are eigenstates of the total Hamiltonian and /3 is the inverse temperature. In 
momentum space it has the dispersion representation [1,2] 

where /(a;) = l/[e^'^ — 1] and the spectral function is 

1 ^ /■ p-^En 



AT 



Because p{k) is real, it is easy to separate the real and imaginary parts of (2): 

p(a;, k) 



ReD,,{k)=r 

Bk^ 

linDii{k) = -7rcoth(-— ) p{k) 



(4) 



To explain what is meant by the Feynman propagator at non-zero temperature it 
is helpful to recall that real-time calculations require doubling the number of degrees of 
freedom [2-9]. Associated with each physical field is an auxiliary field. All propagators 
become 2x2 matrices in the internal space. Thus Dn is one entry in the 2x2 matrix. 
In momentum space it has the form [2-9] 



Dad{k) = Uab 



Dpik) 
-D},{k) 



Ucd (5) 

be 



where Dpik) is the Feynman propagator. The simplicity of this matrix structure indicates 
that the thermal Feynman propagator plays a central role in thermal field theory. To 
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extract a representation for Dpik), use U12 = U21 = [exp(/?|A;°|) — 1] ^^'^ and Un = U22 = 
exp(/3|/cO|/2) U12 to obtain 

^''^^^ exp(/?|/cO|)-l 

When the real and imaginary parts of this relation are compared to (4) the result is 

p{uj,k) 



/oo 
-00 



k^-u (7) 
lmDF{k) = -7re(A;°) p{k) 

Consequently the thermal Feynman propagator satisfies the dispersion relation 

DF{k) = r rfo.— ^^^^^l^^L^ (8) 
^ ^ y_oo k^ - CO + ir]€{k^) ^ ^ 

Using p{—k , k) = —p{k , k) this can also be written 

DF{k)= fc/cj^^^^^^ (9) 
Jo {k^f-u^ + if] ^ ' 

This dispersion relation is not new; it is discussed in [2] for example. It shows that 
the thermal Feynman propagator is an even function of A;^ that is analytic in the region 
Re(A;°)^ > 0. The dispersion relation has exactly the same appearance as at zero tem- 
perature because all dependence on the temperature is contained in the thermal spectral 
function. (In practice it is more complicated than at zero temperature because there are 
no regions of a; where the spectral function vanishes.) 

Since Dp has rather simple properties and plays a central role in the matrix (5), it 
seems worthwhile to ask how Dp can be expressed directly in terms of matrix elements 
of the field operator 4>{x). The answer to this question is given below in (10). The proof 
of (10) is that the Fourier transform produces the defining dispersion relation (8). The 
Appendix contains a discussion of how to extract the operator form for the Feynman proper 
self-energy '^pix). 
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-I3En ^^^^ 



II. OPERATOR FORM FOR Dpix) 
In coordinate space the Feynman propagator is the thermal average of a certain or- 
dered product of two field operators. The ordering is, however, more complicated than the 
usual time-ordering. The result is 

N 

+(iV|0(O) ONi-t) (j>{x)\N) 
The new ordering operation is defined by 

eN{t) = e{t)e{H - e^) - e{-t)e{EN - h) (ii) 

Having the Hamiltonian operator in the argument of the theta function is a compact way 
of representing a projection operator. For example 

e{H-Er,)= (12) 

Ea>En 

Since the ordering depends only on the energy En and not on a particular state \N), it 
could be labelled dE^ix^), but that seems cumbersome. The time derivative of (11) is 
particularly simple: 

j^ONit) = S{t) (13) 

Note that (10) can be written in a variety of forms by inserting a complete set of states 
(12) and regrouping the matrix elements in various ways. In the zero-temperature limit 
(10) reduces to the usual time-ordered product as it should. In this limit the only state 
|A^) that contributes is the vacuum. Since 9{H — E^ac) — 1 and 9{E^ac~ H) — 0, it follows 
that ^Ar(t) 9{t) at zero temperature. 

The remainder of the discussion will prove that the Fourier transform of (10) yields 
the dispersion relation (8) . To do this it is simplest to transform only the time-dependence 
of (10) without changing the space-dependence 

/•oo 

DF{k^,x)^ dt e'''°* Df{x) (14) 



Examine the first term {N\(f){x)9N{t)(f>{0)\N) in (10). Using the time-dependence of the 
Heisenberg field, 0(a;) = exp{iHt)(f){x) exp{—iHt), this can be written 

{N\cl>{x)eMmmN) = (iV|0(f)e^(^--^)* dN{t)m\N) 
the Fourier transform is 

/oo 
dt e*^"*(iV|0(a;)^jv(t)0(O)|iV) = i{N\(j){x) R 0(O)|7V) 
-oo 



where 



„ _ e{H - En) ^ e{EN ~ H) 

= T7^ ^ : r 



+ En - H + ir] + En - H - ir] 
As before, the appearance of the Hamiltonian operator can be replaced by summing over 
a complete set of states. One can write R more compactly as 



oo 



I d-^ ^^^"- 'Hi (15) 
' P - a; + i'qe{k^) ^ ' 



In the same way the Fourier transform of the second term in (10) is 

/oo 
dt e'''°*{N\(l){0)dN{-t)(f>{x)\N) = -i{N\(f>{0) S (f>{x)\N) 
-oo 



where 



oo 



S=l ^ - ) (16) 

' k^ - u) + ir]e{k^) ^ ' 



oo 



The propagator (10) thus has the form 

^ Z k^-uj + irie{k^) 
where 

fN{oJ, x) ={N\cf>ix)5iuj + En- H)cj>iO)\N) 

-{N\(/,{0)d{uj -En + H)(j){x)\N) 
It is easy to see that /jv is the Fourier transform of the commutator: 

df / 

— e-*((iV|0(a')e^(^--^)V(O)|iV) 



-(Ar|0(O)e*(^-^^)V(^)|A^)) 

r°° dt 

J_^^e^-\N\[<Pix),m]\N) 



where the time dependence of was used in the last step. Thus the integrand of (17) 
contains the thermal spectral function: 



The spatial Fourier transform of this is the defining relation (8) and thus proves that the 
operator form (10) is correct. 

To use the operator form systematically one could develop the diagramatic rules for 
perturbation theory. This would require a Wick theorem for the 9n product and an 
extension of the operator approach of Nieves [8] from the time-ordered to the Feynman 
case. It might also be interesting to find the operator form for the finite temperature 
Feynman propagator of quarks and gluons. 
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Consequently 
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APPENDIX A: THE FEYNMAN PROPER SELF-ENERGY 

It is possible to deduce the operator form for the Feynman proper self-energy IIf{x) 
that generates Df{x) . In momentum space the proper self-energy Hpik) is defined by 

The dispersion relation (9) guarantees that Hpik) is analytic in the region Re(A;°)^ > 
and that Imllp{k) is negative when real. Rewrite (Al) in the form 

{-k"^ + m^)DF{k) ^ -l-Upik) Dpik) {A2) 

The coordinate space form of (A2) is the Schwinger-Dyson equation: 

{n + m'')DF{x) = -S\x) - J d^y Uf{x - y)DF{y) {AZ) 

To deduce the operator form for IIf{x) we therefore need to apply (D-l-m^) to the operator 
representaion (10). A single time derivative of (10) gives 



(^4) 



iDF{x) = YX{N\k^) dNit) m\N) 

N 

+{N\m e^i-t) <j>{x)\N)\ 
Because of (13) the terms involving the time derivative of ^iv(^) give 

^ „-/3Epf 

S{t) Y.(N\ [Hx), m] \N) = (A5) 

AT 

which contains the commutator is at space-like separation and thus vanishes by causality. 
The second time derivative is 



AT 



+ {N\4){x)eN{t)(/}mN) (^6) 



+{N\(f>{0) 9N{-t) 4>{^)\N) 



e 



The first term of (A6) gives —i5'^{x) because of the equal time commutation relations. The 
remaing two terms are given by the field equation 

4>{x) = (V^ - m^)(f){x) + J{x) {A7) 

where J = SCj/Scf). Consequently (A6) becomes 

{n + m^)DF{x) = -6\x)-tJ2[{N\J{x) ^N{t) 0(O)|iV) 



AT 

+{N\<p{0)eN{-t)J{x)\N)' 



(AS) 



Z 

This is just the form (A3), li Dp in (A3) were written as a function of x — 2; instead of x 
only, then the integrand of the d'^y integration would heUpix — y)DF{y — z). Comparison 
with (A8) gives 

/ df^ynpix - y)DF{y - z) = iY\{N\J{x) - ^°) 0(^)|iV) 

(^9) 

+{N\<j>iz) 9Miz^ - x^) Jix)\N)\ 

This defines the operator form of the Feynman proper self-energy in direct analogy to the 
definition of the time-ordered proper self-energy [10]. 
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